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ABSTRACT 


The  stress  analysis  of  a  rectangular  tensile  bar  weakened  by  surface 
cracks  Is  formulated  in  tens  of  mapping  functions  of  a  complex  variable. 
Although  an  approximate  solution  of  this  problem  has  been  previously  con¬ 
sidered  by  Tait,  the  present  derivation  specifies  the  proper  load-free 
boundary  conditions  on  the  lateral  edges  of  Idle  rectangle  as  well  as  allow¬ 
ing  for  the  variation  of  both  crack  length  and  length/vldth  ratios  of  the 
plate.  Determination  of  the  stress  function  is  reduced  to  the  solution 
of  a  linear  system  of  simultaneous  equations . 


OSCAR  L.  BOWS 
Mathematician 


APPROVED: 


INTRODUCTION 


In  the  area  of  crack  studies,  a  problem  of  practical  importance  is 
the  simple  geometry  of  a  rectangular  tensile  specimen  weakened  by  surface 
cracks.  In  order  to  apply  such  theories  as  that  of  Irwin1,  it  is  neces¬ 
sary  to  determine  the  stress  distribution  in  the  vicinity  of  the  crack 
root.  Althou^i  the  behavior  of  the  stress  in  such  a  neighborhood  is  known 
qualitatively,  quantitative  results  require  a  formal  solution  of  the  stress 
problem  as  a  whole.  It  is  this  matter  which  is  considered  in  this  report 
for  the  plane  problem  shown  in  Figure  1. 

R.  J.  Taita  has  considered  the  problem  of  determining  the  state  of 
stress  due  to  surface  cracks  in  an  infinite  elastic  strip  of  finite  width. 
Unfortunately,  in  his  analysis,  load- free  conditions  on  the  strip  surfaces 
are  only  partially  prescribed.  In  fact,  a  mixed  boundary  value  problem 
was  assumed  in  his  analysis  which,  in  effect,  is  intrinsic  to  an  infinite 
sheet  containing  a  row  of  periodically  spaced  cracks. 

In  the  present  report,  we  shall  consider  surface  cracks  in  rectangular 
plates  of  finite  length  and  width  with  the  proper  load-free  boundary  con¬ 
ditions  on  the  lateral  edges  of  the  rectangle.  In  addition  to  introducing 
the  desirable  modification  of  Tait's  analysis,  it  is  possible  to  study  the 
additional  parameter,  length/width  ratio  and  its  effect  on  the  stress  con¬ 
centration.  Tbe  polynomial  approximation  method  used  by  the  author3  in 
the  analysis  of  an  infinite  plate  containing  radial  cracks  does  not  appear 
well  suited  to  the  present  problem  due  to  slowness  of  convergence.  Thus, 
an  alternate  approximation  scheme  is  devised. 


CONVENTIONAL  COMPLEX  VARIABLE  FORMULATION 


Muskhelishvlli '  s4  original  complex  variable  statement  of  the  plane 
problem  in  rectangular  coordinates  depends  on  the  representation  of  Airy's 
stress  function,  U(x,y) ,  in  terms  of  the  two  analytic  functions,  cp(Z)  and 
tjf(Z) ,  of  the  complex  variable  Z*=x  +  iy  as 


U(x,y) 


(1) 


With  this  representation,  the  stress  and  displacement  components  can  be 
written  as 


°y  +  ax  "  2  +  cp'^]  *  k  [V(Z)J  *  (2) 

°y  "  ax+  21  V  “  2  [Zq,"(Z)  +  t,(Z)]'  (3) 


3' 


2H  (u+iv)  -  K<f>(Z)-Z  <p'(Z)  - 


(4) 


where  the  pr^me  notation  denotes  differentiation  with  respect  to  Z,  the 
hars  denote  complex  conjugates,  and  U  and  K  are  material  constants.  A  use¬ 
ful  relation  for  handling  stress  boundary  conditions  is  provided,  by 

cp(Z)  +  Z  ?TzT  +  f(Z)  -  i  /(X,  +  i^ )  ds  -  fx(3)  +  if2^  (5) 

o 

where  Xvds  and  Xyds  represent  the  x-  and  y-  components  of  the  force  act¬ 
ing  on  the  element  of  arc  ds  from  the  positive  side  of  the  normal  v . 

Curviliniar  coordinates  can  be  conveniently  introduced  by  considering 
an  auxiliary  complex  Q -plane  along  with  a  conformal  mapping  function 


Z  ■  <u(  0* 


(6) 


The  analytlcity  of  (6)  ensures  the  analytidty  of  the  stress  functions 
Vhen  considered  as  functions  of  the  variable  C  •  New  notation  is  simplified 
by  designating  cp(z)2 cp[uu(C))  as  <p(C),  etc.,  leading  to  such  relationships 
as  <p'(Z)»  <p7(C)/«>  (C)>  etc. 


The  stresses  and  displacements  can  then  be  written  as 

cy  +  ax  -  2  [cp'(c)Mc)  +  vTcTMc)], 


2u  (u  +  iv)  -  k  <Kc)  -  He)  ^TcJ/u’TcT  -  TTcT  • 


(7) 

(8) 
(9) 


When  C  corresponds  to  a  closed  contour  in  the  Z-plane  into  which  the  unit 

circle  C*  a  ■  ei0  is  mapped  by  (6),  the  arc  parameter  S  can  be  considered 
as  a  known  function  of  o,  hence 

tx  (8)  +  i  f2  (S)  -  g  (o).  (10) 

Uie  boundary  condition  on  C  (assuming  the  applied  load  is  given)  can  thus 
be  written  as 


<p(o)  +  <o(o)  p'lerj/w'lcr)  +  Her)'  -  g(o) 


(11) 


OHE  MAFPXJK}  PUKCTIOM 


The  region  occupied  by  the  rectangle  with  surface  cracks  will  be 
denoted  by  3^  with  boundary  T  in  the  the  Z-plane  (Figure  l) . 


The  exact  sapping  function  <u(  0  ie  defined  as  the  mopping  of  the  unit 
circle  and  its  interior  into  the  boundary  t  and  V  respectively. 

By  a  straightforward  consideration  of  the  Schwartz-Chrlstoffel  trans¬ 
formation  ,  it  can  be  seen  that  the  exact  mapping  function  corresponds  to 
an  appropriate  branch  of  the  function 


Z  m  -i  +  l)  j"(l-e2la  C2)(l-e‘21ct  C2)(l-«2iP  C^Cl-*’2^  C2)l  dC. 

L  J  m 
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The  exact  mapping  function  has  eight  branch  points  located  on  the 
unit  circle  necessary  to  the  description  of  the  elgit  right  angle  corners. 
We  define  the  real  parameters  in  (12)  such  that  0<ct  <ps  x/2,  thus,  When 
the  circle  maps  into  the  limiting  case  of  a  rectangle.  The  crack 
roots  are  described  by  the  zeroes  of  cu'(C)  on  the  unit  circle,  i.e.,  C*>  ± 
i.  The  ratios  L/xq  and  Yq/xq  can  be  varied  by  appropriately  adjusting  the 
parameters  a  and  £. 

In  series  form, 

z  »  u<C)  -  -i  JJl  An  C 2x1-1  (13) 


where  the  A&'s  are  real.  The  A^'s  may  be  obtained  numerically  from  a 

straightforward  recursive  formulae  (Appendix  I)  determined  by  comparing 
(12)  and  (13).  The  series  (13)  converges  interior  to  the  unit  circle  and. 
on  any  Interval  of  the  unit  circle  not  Including  a  branch  point. 

Unfortunately,  the  rate  of  convergence  of  (13)  is  slow  on  the  unit 
circle,  thus  leading  to  extremely  tedious  numerical  calculation  in  the 
stress  analysis.  On  the  other  hand,  the  procedure  of  polynomial  approxi¬ 
mation  of  the  mapping  function  utilized  in  Reference  3  is  plagued  by  the 
same  slowness  of  convergence. 


We  adopt  a  different  approach  and  define  for  K21, 

4C,K)«  -i  /(C2+  1)  [(K2- e21aC2) (K2- e-21aC2) (K2- e21V) (K2- e-2iPC2)] 


dC. 

(14) 


Clearly  <o(C,l)  coincides  vlth  (12).  On  the  other  hand,  for  K>1  the  branch 
points  are  moved  out  to  lie  on  a  circle  |  C  |  *  K  >1.  The  cusps  represent¬ 
ing  the  crack  roots  are  still  preserved  on  the  unit  circle  at  +  i.  Thus, 
by  evaluating  (l4)  for  values  of  K  w  1  on  the  unit  circle,  increasingly 
accurate  approximations  T  are  obtained  as  K-*  1.  On  the  other  hand,  the 
convergence  of  the  series  expansions  for  K  >1  will  be  substantially  improved 
compared  with  that  of  (12) .  The  form  of  the  power  aeries  (13)  is  preserved 
in  the  expansion  of  <ju(  C,K) ,  thus, 


»(:.«  - JbA  c2"'1.  («) 

where  the  A  's  are  real  and  are  now  considered  as  functions  of  the  parameter 
K  as  veU  as  a  and  p. 
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DETEH-m&TION  OP  THE  STRESS  FUNCTIONS 


It  is  necessary  to  determine  the  two  stress  functions  <p(t)  end  6(C) 
such  that  both  functions  are  analytic  for  |  C  |  <1  and  the  condition  (11) 
is  satisfied  on  the  unit  circle.  Apart  from  mathematical  elegance,  little 
is  gained  by  any  consideration  other  than  a  direct  power  series  argument 
in  the  present  formulation.  Die  difficulty  due  to  the  zeroes  of  u>*(C)  on 
the  unit  circle  can  be  eliminated  by  writing  (ll)  in  the  form 

<u'(a)  <p(o)  +  <u(c)  <p'(cr)  +  tu'(a)  H®)  ■  ujH®)  6  (<*)•  (l6) 


From  symmetry  considerations, 

*<0  -  1  1  nil  *ta  f*'1’  <1T> 

♦(C)  -  i  1  J^n  C2”'1,  (18) 

where  -the  coefficients  a  and  b  are  real.  If  the  stress  state  were  that 

of  uniform  tension,  l.e. ,  no  surface  cracks  were  present,  -then  <p(C)  ■  T  X 
u)(C)/4  8X1(4  HO  “  T  «i(C)/2.  Thus,  we  assume 


9(C)  -  T  HOA  +  Vj.  (0, 

t(C)  -  ®  HC)/2  +  (0, 

Where 

Clearly, 

<p1  (a)  +  ufa)  9{(<t)  +  ("'(a)  J^aJ  - 


i^Cc)  - 1  *  JL  °n  c2”"1* 

♦i(0  •  1  1  A  Pu  C2"’1. 


(19) 

(20) 

(21) 

(22) 


Substitution  of  (21)  and  (22)  into  (23)  and  comparing  coefficients  of  equal 
powers  of  o  leads  to  the  two  sets  of  conditions 


Jl  A  (2n-l)  a 


21  a  (2n-l)  A 


n?p+2  An  (2n"1)  an-p-l  +  n*p*-2 


an  Vp-1 


‘  2l  An  <2n‘1)  V-n  '  ‘  V  p  '  °*1*2' -  (25) 

Hie  set  of  equations  25  provides  for  a  direct  determination  of  the  p  *s 
once  the  a  's  have  been  found  from  the  linear  system  of  simultaneous nequa- 
tlons  24.  n 

Hie  applied  load  will  he  expanded  on  the  basis  of  the  exact  geometry 
and  loading.  Hius,  if 


8  (a)  -  |  [w(<0  +  -  jJC"  (^2K-1 

then  it  can  easily  be  shown  (Appendix  2; ,  that 

°ac-i  -  -  ™  Jl  ao  ■l“  t20-1)  P 


(26) 


+  ~  J*  '  An  sin  (2n-l)  0  sin  (2K-1)  OdO,  K-<*±3,±2, - (27) 

P 


«/2 


(28) 


For  conputatlonal  purposes  an  alternate  form  of  C2K_1  is  now  presented. 

c2K-i  -  -  ^  [*te)  -  x(®>]  8in  (sk-1)  ®  *  0, 

where  X(o)  is  the  real  part  of  cu(e 10 ^ .  Integration  by  parts  yields 

2iT  f  ^  <2%  ___  a  4  a  (29) 


c  .  . ari  .  r 

°2K-1  (2K-1JX  ^ 

Using  (12)  to  confute  dX/dO,  it  follows  that 


cos  (2K-1)  OdO. 


-8. 


n  m  .  r _ Mg  9  .SO ft  lags!}  _ 

-1  (£k-1)k  £  (cos  2a  -  cos  2©)^' 2( cob  2£  -  cos  2©)1/2 

From  the  definition  of  the  d^'s  in  (23),  it  follows  that 


*/2 


.sag.,  ft  -sag  (3fcl)ft,..&.,ft.. 


(30) 


S  “  1  a&.  An  C2(!*n)-3'  p  -  O'*1' *2, -  (3D 


When  approximate  mappings  13  are  used,  the  loading  conditions  must 
be  considered  as  approximated.  On  the  other  hand,  as  K-+1  both  the  geome¬ 
try  and  the  applied  load  converge  to  the  exact  problem. 


DETEIMafcTIQN  OP  THE  STRESS  CXJNCEinmTION  IK  THE 
VICINITY  OF  THE  CHUCK  HOOTS 


The  stress  concentration  in  the  vicinity  of  either  of  the  crack  roots 
is  the  prime  aim  of  the  analysis.  If  polar  coordinates  (r,y)  are  Introduced 
with  the  origin  chosen  as  the  crack  root,  the  stresses  vary  to  the  first 
order*  as 


a  «  L(7)/YF*.  (32) 

The  determination  of  L (7)  in  terms  of  the  present  analysis  will  now  be  con¬ 
sidered. 

Assume  that  C  is  the  point  on  the  unit  circle  which  maps  into  the 
crack  root  Z  .  Then  u)'(  $  has  a  simple  zero  at  C«  C  •  Furthermore,  00(C) 
la  analytic  £n  the  neighborhood  of  Q  on  the  Rlemann  surface  defining  the 
mapping  function,  thus, 

JL  Sn<C-to>n-  (33) 

By  reversion  of  series, 

£-  !o  -  <z-z0>n/2’  <»> 


*  L(7)  is  essentially  proportional  to  Irvin's  factor  K. 


,9- 


where 


2 

Sg  ■  1,  etc. 


thus, 


c-c»“#<z-2°>1/2- 


(35) 


If  ve  introduce  polar  coordinates. 


Z-Z  =  r  ei7. 
o  9 


then 


(.r,£  .ir/2  .  ^  .W/2 


V% 


(36) 


(37) 


In  order  to  utilize  these  approximations  to  examine  Idle  stress  field 
local  to  the  notch  root.  It  is  clear  from  (7)  and  (8)  that  the  character 
of  cp(C)  and  $(C)  in  the  neighborhood  of  Is  required.  If  ve  define 


where 


*(C)  -  -  ?(i/c)  -  «(i/c)  9'(C)/«'(C) 


f(i/o  = 


(38) 


(39) 


and  substitute  into  (ll),  it  can  be  shown  that  <p(C)  is  analytic  in  the 
neighborhood  of  C=CQ  whereas  $(C)  has  a  single  pole  in  this  neighborhood. 


It  follows  fairly  directly  from  (7)  that 


ox  +  Oy  «  KL  \ 


-17/2 


V37Eo5  V7 

b 

From  (8),  it  follows  from  a  somewhat  more  difficult  calculation,  that 

-/2“e-i7/2 


(«K» 


«  -  o_  + 


w 


9'(c0)  ^•(Co)  ^21y 
T2 - 5^ — e 


.  (41) 
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In  the  calculation  of  (40)  and  (4l),  the  key  quantity  is  cp'CCq). 
Whereas  the  calculation  of  u)"(  r  )  can  be  made  directly  from  the  closed  form 
of  the  mapping  function,  the  determination  of  9'(C  )  requires  the  solution 
of  the  system  of  equations  24. 


SUMMARY 


Die  stress  analysis  for  a  rectangular  tensile  bar  weakened  by  surface 
cracks  has  been  presented.  A  mapping  function  which  permits  variation  of 
the  crack  depth  and  length/width  ratio  of  the  plate  Is  utilized.  Determina¬ 
tion  of  the  stress  field  in  the  neighborhood  of  the  crack  root  requires  the 
calculation  of  <p'(C)  vhere  C  is  the  point  on  the  unit  circle  correspond¬ 
ing  to  the  crack  root.  A  procedure  for  systematically  approximating  cp'(c) 
is  suggested  by  the  limiting  process  as  K-frl  in  the  mapping  function. 

Calculations  are  underway  for  pertinent  values  of  the  parameters  and 
the  results  will  be  presented  in  a  subsequent  report. 
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APPENDIX  I 


RECURSIVE  70MJIAE  TOR  THE  MAPPING  COEFFICIENTS 


Fran  (l4). 


co'(C,K)  -  -i  (C2+1)/(C1;-2K2C2  cos  aa+K1*)1^2^21^  008  2^+K^) 


1/2 


Let 

tf(C,K)  -  -i  ^  C211"2. 

Then, 

^K2  (-2  cos  2{3  »2  cos  2a) 

+  C4  (2K4+te4  cos  2a  cos  2p)  +  r2  (-2K^  cos  2a  -  ac6  cos  2p) 

♦»?}  {& • 

Let 

2e^  -  -2  cos  2a  -  2  cos  20 
2cg  »  2+4  cos  2a  cos  20 

a,  -  m® 

C3  ■  2Mi  Mg 

C4  -  sm^  +  Mj,2 

C5  -  2M^  Mj^  +  2Mg  etc. 

Then, 

(C4  +  2  ^+1)  -  ^C®+2  K2^  (^  +  2  K4  *2  C4  +  2  K6  ^  C2  +  K8 J 


Equating  coefficients  of  equal  powers  of  C  yields 
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K8  -  1 

K8Cu  +  2K6c1C3  +  2KU€2C^-1 
K8  C5  +  2  K6  €1C4+2K*V€2C3  +  2K2e1C2-0 

K8Cg  +  2K8c1C5  +  2K1*C2C^+2K?c1C3+C2-0 

etc. 

Thus, 

^  -  -l/K4 
Mg  -  -K4  C3/2 

^  -  -K4[c4/2  -M22/2j 
M4  "  -k4[c5/2  -  ^  Mjj 
etc. 


Thus,  in  the  expansion 


»  (C,K)  -  -  i 


An  -  y(2n-l). 
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APPENDIX  n 


EXPANSION  OF  THE  LQADQB  JUNCTION 


It  vill  be  recalled  that 


S 


sM  •  i  J  (X,  +  i  *,)  33 


We  define 


g 


(a)  .a,  [jc0  *  X  *  0,  y  -  jrJ] 
.-St0,[x  -  Jt0,  O*y*y0] 


etc. 


!Qius, 


g  (a)  -  |  j»  (a)  +  u>  (ajj  •  0,  0  <  9  *  p 

-  -  aXQ  -  |  |w  (a)  +  aTTa^  ,  p  s  Q  *  x  -p 


0,  x-psOsJx  +  p, 


-  TX©  -  |  j»  (<0  +  “*?)J  > 


-0,  2x-px©s2x 


+  x  s  »  s  2 


x  -pj  , 


II  we  set 


g  (or)  - 1  [.  (ff)  ♦  in  lo;j 


% ' K- 


then 


14- 


<r‘  b  l  "*  ^-xo '  i  [“M  ♦  ^]J 

|x0  -  5  ♦  *.)]j 


-1X0 

e  d0 


46. 


R+0 


It  follows  from  syametiy  that 


»  0,  K  even 
T  r*“^  I 

<*  -  I  l 


P 


Again  from  symaetry. 


-*0  •  5  ^ct)  +  «r»yj 


,4jro 

e  d®,  K  «  ±  1,  ±  3#— 


Bin  KQ  dO,  K  «  ±  1,±  3, - 


Finally, 


<W  *  O'  K  -  0,  *  1,  ±  2,  - 


<k-l-*¥  C°aATl)  P  Jl  An  ,ln  &-1'  ^ 


21T  */2  * 


+  *  !*  n&.  An  #in  (2np>1)  ®  8in  (2K*1)  ®ao*  K-°,  *  1,*  2,- 
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